Win?. 


CLEMSON  UN IV  S  C  DEPT  OF  MATHEMATICAL  SCIENCES  P/f 

THE  INITIAL  BUST  CYCLE  OF  A  DISCRETE-TIME  81/6/1  8UEUE  IN  WHICH  ■CTCl. 
JUN  79  D  L  MINH  N00018-75-C-04S1 

UNCLASSIFIED  NUO  NL 


THE  INITIAL  BUSY  CYCLE  OF  A  DISCRETE¬ 
TIME  GI/G/1  QUEUE  IN  WHICH  THE 
INTER-ARRIVAL  TIMES  ARE  NOT 
NECESSARILY  IDENTICALLY  DISTRIBUTED  , 

Do  „e  Minh* 


Clemson  University 

Technical  Report  #319 
JunaSP.  1979 


*On  leave  from  The  New  South  Wales  Institute  of  Technology,  P.0.  Sox  123, 
Broadway,  NSW  2007,  Australia. 


Research  supported  in  part  by  THE  OFFICE  OF  NAVAL  RESEARCH  Contract 
N0001 4-75-C-0451 j 


DT! 


L  • 

At 

A  1  A  -V  ,  ' 

J\  Cj 

:a r  a  , 

Apptov«d  ‘->i  i  •»>  ; 

B 

Distribulion  Urdu);.:  d  • 

THE  INITIAL  3USY  CYCLE  OF  A  DISCRETE-TIME  GI/G/1  QUEUE  IN  WHICH 
THE  INTER-ARRIVAL  TIMES  ARE  NOT  NECESSARILY  IDENTICALLY  DISTRIBUTED. 


Do  Le  Minn , 

Department  of  Mathematical  Sciences, 


Clemson  University, 

Clemson,  SC  29631,  U.S.A.  (’ 


ABSTRACT . 

This  paper  studies  a  generalization  of  the  discrete-time 
GI/G/1  queueing  system.  Here,  the  inter-arrival  times  are  nor 
necessarily  identically  distributed.  A  recursive  scheme  is 
derived  to  obtain  the  joint  distribution  of  the  duration  of 
the  initial  busy  period,  the  duration  of  the  first  idle 
period  and  the  number  of  customers  served  during  the  initial 
busv  oeriod. 
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Humber  of  customers  served  during  the  initial  busv  ceriod. 
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The  initial  busy  cycle  of  a  queueing  system  has  always  been 
a  subject  of  interest.  It  is  of  great  importance  in  the  study  of 
systems  in  which  the  cost  associated  with  unproductive  idle  period 
is  high  and  thus  the  server  has  to  be  'switch  off'  or  transferee 
to  another  system  whenever  the  original  system,  is  empty  for  the 
first  time.  Most  of  the  results  obtained  in  the  literature  for 
the  initial  busy  cycle,  however,  are  not  very  suitable  for  computer 
applications.  For  the  GI/G/1  queue,  they  are  presented  in  terms 
of  Laplace-Stiel jes  transforms  involving  complicated  operators 
(  [81/  ecu.  5.13).  For  the  M/G/l  queue,  they  contain  the  roots  of 
a  certain  functional  equation  within  an  unit  circle  {[3],  ecu.  25; 
Obtaining  numerical  information  from  these  results  is  not  an  easy 
task  and  has  never  been  discussed  thoroughly  in  the  literature. 

In  this  paper,  we  study  a  discrete-time  single-server  queueing 
system  in  which  the  inter-arrival  times  are  mutually  independent 
but  not  necessarily  identically  distributed;  the  service  times 
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served  during  the  initial  busy  period.  This  scheme  is  distribution 
f ree ;  that  is,  it  requires  no  specific  distributional  assumptions 
for  the  random  variables  underlying  the  model.  This  scheme  _s  suite 
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for  systems  having  a  very  general  initial  condition.  When  the  system 
is  initiated  by  the  arrival  of  a  customer  to  an  empty  queue,  this 
scheme  will  enable  us  to  obtain  information  about  the  normal  busy 
cycle.  For  the  queueing  system  in  which  there  is  a  random  'sec-up' 
time  for  customers  who  arrive  when  the  server  is  idle ( [7] , [3] )  this 
scheme  is  also  applicable  as  the  sec-up  times  only  affecc  che 
initiation  of  the  busy  cycle. 

For  the  actual  waiting  time  of  each  customer  in  this  model, 
another  recursive  scheme  has  been  reported  in  [9] . 

Let  us  restricc  ourselves  to  the  discrece-time  system  observed 
at  equally  spaced  epochs  0,  1,  2,  n,  ...  All  events  of  the 

system  namely  arrivals,  transfers  from  queue  to  service  and 
departures  are  assumed  to  occur  at  instants  immediately  prior  to 
these  epochs.  We  assume  that  initially  the  system  is  not  empty 
and  that  there  are  (kQ>i)  customers  in  it  at  che  epoch  r.=0 
and  the  (kg+i}th,  0<g-2)th,  ...,  kth,  ...  customers  arrive  at 
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where  0<~. 


<c,<...  We  also  assume  that  the  inter-arrival  times 
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identicallv  distributed.  We  write 
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the  service  time  of  each  customer  be  3.  We  write 
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The  customers  are  served  in  order  of  their  arrivals  arc 


no  limit  on  the  waiting  room. 
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First,  let  V  (n>0)  be  the  virtual  waiting  time  at  epoch  n. 

~n 

Let  A  (niO)  be  the  time  difference  between  n  and  the  time  of  the 

first  arrival  in  the  interval  (n,»).  Let  C  (n>0)  be  the  total 

~n 

number  of  customers  being  in  the  system  at  the  epoch  n=0  or  arriving 

at  the  system  during  the  interval  (0,n].  We  shall  study  the 

multi-variate  Markov  chain  {V  ,A  ;C  }  . 

~n  ~n  ~n 

It  is  well-known  that  any  stochastic  process  can  be  character¬ 
ized  as  a  Markov  process  if  the  full  set  of  random  variables  needed 
to  specify  the  state  of  the  process  is  employed.  For  the  queueing 
system  studied  in  this  paper,  the  instantaneous  state  of  the  process 
is  completely  characterized  by  the  bi-variate  {Vn,A  }  .  In  this 
paper,  however,  instead  of  studying  the  bi-variate  {V  A  }  ,  we 
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study  the  tri-variate  {V  A  ;C  }  .  We  shall  show  that  this  will 
-1  ~n,  -n  ~n 

enable  us  to  resolve  our  problem  more  fully  (See  also  [6]  ,  [7]  ,  [9]  )  . 

Now  let 
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Pr{vn=j , An=l;Cn=k}  (n>0, j>l, i>l,k>l) 
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Let  us  now  consider  a  second  queueing  system,  having  exactly 

the  same  description  as  the  one  above  but  admitting  no  more  customers 

whenever  the  system  is  empty  for  the  first  time.  It  is  equivalent 

to  the  initial  busy  period  of  the  original  queueing  system.  The 

second  system  will  have  the  same  Chapman-Xolmogorov  difference 

equations  as  (4) -(6)  excect  that  the  term  v*  .  ,  ,s.a'.  does  not 

0 , 1 ;  .<  -  j.  j  L 

appear  in  ( 6 )  . 

Let  the  notation  for  the  second  queueing  system  be  those  for 
the  first  system  together  with  the  additional  symbol  '  .  We  say 

that  the  chain  { ;  has  {0,i;k  <->1 , k >1 }  as  absorbing  stares. 
Cox  and  Miller  [2, p.236]  suggest  this  approach  as  one  suitable  men bed 
to  investigate  the  distribution  of  the  busy  period  of  z'r.e  M/G/l 
model  but  do  not  elaborate  any  further.  (See  also  [4]  and  [5]). 

From  (1),  (2),  (5) -(6),  with  the  term  v”  .  ,  omitted, 

after  some  routine  calculations,  we  obtain  the  following  relation 
between  the  generating  functions 
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do  so.  It  is  therefore  equal  to  the  joint  probability  that 
the  initial  busy  period  of  the  original  system  ends  at 
epoch  (n+1) ,  there  are  k  customers  served  during  this  period 
and  the  duration  of  the  first  idle  period  is  i^-l).  If  we 
denote  this  by  b1?*^  ,  ,  we  shall  have 
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This  relationship  indicates  that  the  busy  period  studied 
in  this  paper  is  what  Cohen  [1,  p.  284]  calls  the  strong  busy 
period.  A  strong  busy  period  continues  when  an  arrival  occurs 
at  the  instant  of  departure  of  the  last  customer  in  the  system. 
In  other  words,  two  strong  busy  periods  are  separated  with 
probability  one  by  an  idle  period  of  non-zero  duration. 

From  (7) ,  taking  the  generating  function  with  respect  to 
time  yields 
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If  0<U!*|t|4l,  k^l,  letting  z  =  i/t  in  the  above  equation 
we  obtain 


(10) 
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From  (10)  ,  denoting  by  the  coefficient  of  zu  in  ^(z) 

and  by  a™'n  the  coefficient  of  zu  in  n^_m-lv(z)  ,  by  induction, 
we  get 
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Equating  the  pc  it.  tire-degree  coefficients  of  t  in  the  above 
equation  yields 
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Further,  equating  the  coefficients  of  •  yields  the 
following  relations  which  will  allow  us  to  calculate  easily 
the  values  of  b1?  ,  (n>l,  1,  k^2)  from  the  uwe,:  initial 
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